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A method is described for the reconstruction of the amplitude and phase of the object exit wave 
function by phase-plate transmission electron microscopy. The proposed method can be consid- 
ered as in-line holography and requires three images, taken with different phase shifts between 
undiffracted and diffracted electrons induced by a suitable phase-shifting device. The proposed 
method is applicable for arbitrary object exit wave functions and non-linear image formation. Ver- 
ification of the method is performed for examples of a simulated crystalline object wave function 
and a wave function acquired with off-axis holography. The impact of noise on the reconstruction 
of the wave function is investigated. 



INTRODUCTION 

Recent advances in the development of physical phase 
plates have opened new imaging capabilities in transmis- 
sion electron microscopy (TEM). There are different ap- 
proaches to realize such phase plates, for example the 
electrostatic Boersch phase plate [U [2] , drift-tube phase 
plate [2], or carbon-film-based Zernike or Hilbert phase 
plates [H [S]. Except for the Hilbert-type phase plate, 
all these phase plates have in common that they aim 
at shifting the relative phase of the diffracted and un- 
diffracted electrons. In the weak-phase object approx- 
imation phase plates provide contrast enhancement at 
small and intermediate spatial frequencies compared to 
normal bright-field images. However, in any TEM image 
only the absolute square of the image wave function is 
observed with the consequence that all phase informa- 
tion is lost. Moreover, due to imperfect lenses the image 
wave function is modified by aberrations compared to the 
object wave function, meaning that information at spa- 
tial frequency u of the object is transferred with a differ- 
ent amplitude and different phase into the image. This 
also occurs in phase-plate imaging, even in combination 
with an aberration corrector. A constant phase-contrast 
transfer function of one would result in a constant am- 
plitude contrast transfer function of zero. To retrieve 
the full phase and amplitude information of the origi- 
nal complex object wave function is therefore a general 
problem in transmission electron microscopy. Since Ga- 
bor's [SI [7] proposal for holography being used in TEM, 
several methods were proposed and several were realized 
experimentally. For material science specimens off-axis 
holography is a well known and quite useful method for 
amplitude and phase reconstruction j8] [9] . However, the 
small area which can be reconstructed is a disadvantage 
of off-axis holography. This restriction can be alleviated 



by a Lorentz lens but at the expense of spatial resolution. 
Another technique relies on through-focus series, where 
several images at different defoci need to be taken pTOT 
fl"2] . This in-line holography method requires exact de- 
termination of aberrations and defocus of each image for 
correct reconstruction. The method proposed by Danev 
and Nagayama [13] uses a Zernike- or Hilbert-type phase 
plate to reconstruct the object wave function for weak- 
phase objects, which are commonly investigated in biol- 
ogy. Images taken with phase shifts induced by a phase 
plate we refer to as 'phase-contrast images' in the fol- 
lowing. According to Danev and Nagayama, phase re- 
trieval of weak-phase objects requires a phase contrast 
and conventional TEM image. The image wave function 
of a weak-phase object is obtained by adding the Fourier- 
transformed conventional image and complex conjugate 
of the Fourier-transformed phase-contrast image. The 
object wave function is obtained by subsequent correc- 
tion of the image wave function with the wave transfer 
function. 

In this report we present a method for the reconstruc- 
tion of the complex object wave function of arbitrary 
objects which do not have to be necessarily weak-phase 
objects. The method requires the use of a physical phase 
plate or any other device which is capable of imposing 
at least two different relative phase shifts between un- 
diffracted and diffracted electrons. Three images need 
to be taken at three different arbitrary phase shifts, for 
example at — 7r/2, and at ir/2. All other parameters 
like defocus, Cg-value and exposure time must be kept 
constant in all three images. The technique allows a full 
analytical reconstruction of the object exit wave function. 
The complex wave function can be reconstructed without 
any restrictions because the reconstruction does not only 
work for weak-phase, but also for strong- phase and/or 
strong-amplitude objects. It also does not matter if the 
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images are dominated by linear image formation or if 
they contain non-linear contributions. The validity of the 
approach is demonstrated by an experimental wave func- 
tion, which was obtained by off-axis holography. More- 
over, non-linear image formation was taken into account 
to calculate images of a simulated crystalline object wave 
function. The images were calculated for different phase 
shifts of the undiffracted electrons. The reconstruction 
process based on such three images yields the local am- 
plitude and phase which coincide well with the original 
wave function. 



EXPERIMENTAL TECHNIQUES 

To obtain a real test wave function for the recon- 
struction process, off-axis holography in a transmission 
electron microscope was applied to reconstruct the lo- 
cal amplitude and phase of a test sample consisting of 
nanoscaled Pt-particles dispersed on a thin amorphous 
carbon (a-C) film. For off-axis electron holography, a 
Philips CM200 FEG/ST was used which is equipped with 
a Mllenstedt biprism, i.e. a thin wire, in the selected- 
area aperture holder. The electrostatic potential at the 
biprism was set to approximately 180 V corresponding to 
an interference fringe distance of 0.2 nm. Holograms were 
obtained by inserting the sample into the beam in such 
a way that it is positioned on one side of the biprism 
with the specimen edge approximately parallel to the 
biprism. The object wave is transmitted through the 
sample whereas the reference wave propagates through 
vacuum in the hole of the TEM specimen on the other 
side of the biprism. Kinematical diffraction conditions 
prevail for the test sample. The TEM image simulation 
routine of the STEMSIM program [H] written in Matlab 
(The Mathworks Inc.) was extended with the possibility 
to define an additional phase shift 4>pp between diffracted 
and undiffracted electrons. It was also complemented by 
a software package to implement the off-axis holography 
reconstruction as well as the phase-plate reconstruction. 
More experimental details and the off-axis holography 
reconstruction procedure are outlined by Lehmann and 
Lichte [15] , 

WAVE FUNCTION RECONSTRUCTION BY 
PHASE-CONTRAST TRANSMISSION 
ELECTRON MICROSCOPY 

The information of the object after transmis- 
sion through the sample is encoded in the object 
wave function given by ^(r) = \^Jq\ exp(i(f>o(r)) — 
\^d{t) \ exp(i4>z)(r)) which is composed of the local am- 
pltiudes and phases of the unscattered \&o arj d the 
diffracted electron wave ^d- The object wave function 
is modified by the properties of the imaging lens system 



which imposes a phase shift x( u ) = n{Z\u 2 + ^CsX 3 u 4 ) 
on the object wave. The wave aberration function \ de- 
pends on the spatial frequency u and on the objective 
lens aberrations such as defocus Z and spherical aberra- 
tion coefficient Cg. With a phase plate in the back focal 
plane of the objective lens we can change the phase shift 
of the undiffracted electrons </>o(^0 at u =0 in the case of 
a Boersch phase-plate by varying electrode voltage. To 
consider the effect of an ideal phase-shifting device on im- 
age formation, the relative phase shift between diffracted 
and undiffracted electrons is taken into account by in- 
cluding 4>pp in the wave aberration function. 

ti^O : X {u) = tt{Z\u 2 + I/2C S AV) _ | ^ 

u = :x(u) = (f) P p 

Image formation for arbitrary strong scattering objects is 
commonly described by the transmission cross-coefficient 
(TCC) formalism [16]. The intensity (of reflections) in 
the Fourier-transformed image at a spatial frequency u 
is given by 

I(u) = J $(u + u')$* («') T(u + u', u')du' (2) 

with the wave function in the diffraction plane (u). The 
transmission cross-coefficient T(u+u',u') given by 

T(u+u',u') = T{u",u') = E(u" ,u'yexp[i-(x(u")-x(u'))] 

(3) 

which takes into account the phase shift induced by the 
objective lens aberrations x(u). The effects of partial 
spatial and temporal coherence arc described by envelope 
functions contained in E. More details on the explicit 
form can be found in |16[ 117] . 

The expression given by Eq. ^ can be separated into 
the linear contributions (i.e. the interference between 
undiffracted and diffracted electrons u'=0, u'=-u) and 
non-linear contributions (i.e. the interference between 
diffracted electrons (v! 7^' 0, v! 7^ — u). 

I(u) = *(m)** (0) T(u, 0) + (-u)*(0)T* (-u, 0) 

+ / V(u + u')y*(u') T{u + v! ,u')du' + I'(u) 

J u' —u.u'O 

(4) 

I'(u) represents the sum of all contributions to the image 
intensity which are incoherent with respect to the un- 
scattered electron and are not considered by the integral 
in Eq. Q. We will comment on I'(u) below. The inte- 
gral in Eq. Q can be neglected for weak-phase objects 
which yields only negligible contributions to the image 
intensity. Non-linear contributions in images of strong 
objects prevent the application of the method presented 
by Danev and Nagayama [T3] for wave function recon- 
struction which is extended in this work for arbitrary ob- 
jects. Our method outlined in the following subchapters 



relies on the obvious fact that a phase shift imposed on 
the phase of the unscattered electrons by a phase plate 
changes only the intensities of the linear contributions 
leaving the non- linear ones unaffected. I'(u) represents 
contributions to the image intensity which are incoherent 
with respect to the unscattered electrons. This applies to 
all inelastically scattered electrons because a (stationary) 
phase relationship does not exist between unscattered 
and inelastically scattered electrons. We note that inelas- 
tically scattered electrons may produce interference pat- 
terns as shown by energy-filtered off-axis holography |18j . 
However, these interference patterns cannot be exploited 
for wave function reconstruction by phase-plate imaging 
because zero-loss and energy-loss interference patterns 
add up incoherently. 

In the following two sections we describe reconstruc- 
tion formalisms for object wave function in real and re- 
ciprocal space. To give an overview over the steps needed 
for reconstruction an illustration is given in Figure [T] 
Starting with the three acquired images with different 
phase shifts applied by the phase-plate, the two differ- 
ences are calculated after an alignment procedure to cor- 
rect for specimen drift. From these differences the image 
wave function (IWF) is calculated and corrected with the 
wave transfer function. The result is the complex form 
of the object exit wave function. 



Reconstruction in real space 

The reconstruction of the image wave function requires 
three images taken with the same imaging parameters, 
but different phase shifts of the undiffractcd beam. If 
not only the image wave function but also the object 
wave function needs to be retrieved, the aberrations of 
the imaging system need to be known. A wave trans- 
fer function correction can be carried out as described 
e.g. by Lehmann and Lichte [TS]. If a Cs-Corrector 
is available Cs and defocus should both be set to the 
smallest possible value for all three images, which might 
make a wave transfer correction unnecessary depending 
on the slope of the wave aberration function in the range 
of spatial frequencies used for reconstruction. The im- 
age intensity is determined by the sum of linear and non 
linear contributions. Neglecting aberrations for the mo- 
ment, as they can be corrected after the reconstruction, 
we can describe the image intensity in real space by the 
following equation 



i*o-*D(*or 



(5) 



contains the phase and amplitude of the un- 
diffractcd electrons and contains the contributions 
of all diffracted electrons at location r. In the three re- 
quired images, 4>pp for example may be set to 0, it/2 
and — 7r/2. Although the phase shifts can be arbitrary, 
these particular phase shifts were chosen to simplify the 
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FIG. 1. Schematic illustration of the required reconstruction 
steps for the presented method. Input parameters are three 
aligned images and the complex wave transfer function. Image 
wave function (IWF) is calculated and corrected by the WTF. 
Amplitude and phase of the exit object wave function are the 
result. 



following equations. Moreover, the signal-to-noise ratio 
will be better if the difference between minimum and 
maximum phase shift is ir because we see contrast rever- 
sal between the two extremes. Note, that the contrast 
enhancing properties of phase-contrast imaging with a 
physical phase plate, rely totally on weak-phase objects 
and the associated approximation. For strong-phase ob- 
jects, the contrast can decrease for certain <f)pp. It is 
important for the image acquisition that the exposure 
time is identical for all three images, and the require- 
ments for experimental conditions apply as discussed in 
chapter 5. In the next step the difference between each 
of the two phase-contrast images and the " conventional" 
image is calculated. The phase of the undiffracted part of 
the wave V^o is set to zero, because only phase differences 
are of interest, while absolute phases are not observable. 
The result is the difference between the linear interfer- 
ence terms only, as all terms which are independent of 
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the phase shift of the undiffracted beam are eliminated. 
The following equations describe these intensity differ- 
ences 

vf(r)=|* |-|*D(r)|e^ 

*<=br/2)(0 = |*o| e (± " /2) - e #(f) 

£(0^/2) (0 = *(r)**(0 - *(±*/y(f)Vl±*/2)(' f ) 

= |* | 1*0(01 {e im {l ±i) + e(-^»(l T 0) 

(6) 

with the amplitude of the undiffracted wave |*o| and the 
amplitudes and relative phases of the diffracted part of 
the wave function |*£>(r)| and <fi(r). It should be noted 
that 4>{r) denotes a phase linked with the exit wave func- 
tion, whereas the phase shift induced by the physical 
phase-plate is always denoted as 4>pp. Since the phase 
shift tfipp was selected to be ir/2 and — 7r/2, its effect is 
contained in the signs and complex number i in the brack- 
ets of the equations. The equations contain linear image 
contributions because all non-linear contributions as well 
as incoherent (with respect to the undiffracted electrons) 
contributions are not affected by the phase shift that was 
imposed on the undiffracted wave. Constant terms which 
are independent of the phase shift of the undiffracted 
electrons are eliminated. The system of equations can be 
solved, if a third equation is available which is linearly 
independent and contains separately \^q\ and |*£>(f)| ■ 
So far we can solve the system of equations: 

-(-D(0,-tt/2) + «-D(0,+7r/2)) = 

e l<t> (l -i + i-l) + e -i *(l + i + =2e -i *(l + i) 
I 

-(*-D(0,-7r/2) + -D(0,+7r/2)) — 

e**(i + 1 + 1 + i) + e -i *(i - 1 + 1 - i) =2e i0 (l + i) 
with 7 = |*o I |*o| 

(7) 

for 4>{r) and for |\&£)(r)| dependent on |*o|, which are 
given as: 

= I D(Q-k/2) +iD(0,+ir/2) 
]] iD(p,—. ff /2) + -D(0,+7r/2) 

4=±U n ( D < Q >-*™ +iD < > + "V) (8) 

2 «-D(0,-7r/2) + - D (0.+7r/2) 

l * Dl = 2(lli)\*o\ {iD{0 --" /2)+D( °- + " /2)) 

The phase <fi(r) of the diffracted electrons is indepen- 
dent of | 1 1 but it cannot be translated directly into the 
phase of the exit wave function, because it has the form 

*(r) = |*o|- hMOl ( 9 ) 



Different approaches are possible to obtain a third 
equation for | *o | , but they all require some knowledge of 
what is seen in the image. In the following we present a 
method which requires a small vacuum region in the im- 
age where the square of the exit wave function is unity. 
We can then normalize the reconstructed image wave 
function for an illumination with a plane wave with unity 
amplitude by setting |^o| = 1 and dividing the whole 
image set by the vacuum intensity. If a vacuum region 
is not present in the image due to experimental restric- 
tions, it may be possible to find other boundary condi- 
tions for the system of equations. For example if the 
amplitude in a certain area (i.e. carbon film or another 
homogenous area) is known. It should be noted that the 
vacuum could be captured by a fourth image away from 
the region of interest, but the dynamic range and/or bit- 
depth (1 2-bit, 14-bit etc.) of the CCD-camera might pre- 
vent such an image at exactly the same conditions due 
to over-exposure. If the reconstruction is carried out in 
real space the full spectrum of spatial frequencies con- 
tained in the image is used for reconstruction. For some 
applications, especially for images taken without a C5- 
corrector, it might be useful to reconstruct only up to a 
certain spatial frequency u. It should be noted, that the 
same reconstruction method can be carried out in Fourier 
space, with the exception of u=0, for which |* | has to 
be inserted. During the reconstruction process the spa- 
tial frequency interval can be easily limited by use of an 
aperture function. After successful reconstruction of the 
image wave function (IWF), the object wave function can 
be obtained by taking the wave transfer function (WTF) 
into account in analogy to off-axis holography [151 where 
the Fourier-transformed image wave function is divided 
by the wave transfer function in Fourier space. Since 
aberrations should be constant during the acquisition of 
the three images, only one set of parameters needs to be 
determined. It is likely that this can be done with im- 
proved precision from three images, if of course the phase 
shift is taken into account. If an aberration-corrected mi- 
croscope is used, it is possible to correct all aberrations 
as precisely as possible. Since phase contrast is produced 
by the phase plates, it is not necessary to have other than 
residual aberrations for phase contrast transfer |19j . In 
this case no WTF-correction is needed if the wave aber- 
ration function is flat over the reconstructed part of the 
Fourier space. 

Reconstruction in Fourier space 

Based on the same principle the entire reconstruc- 
tion of the object wave function can also be carried out 
in Fourier space. Again three images of the specimen, 
taken at the same imaging parameters but at different 
phase shifts of the undiffracted electrons are required. 
Although these phase shifts can be chosen arbitrarily, we 
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once again choose 0, tt/2 and — tt/2 to simplify the follow- 
ing equations. Unlike the real-space method presented in 
the previous chapter, the reconstruction in Fourier space 
does not rely on a vacuum region in the images of the 
specimen themselves. It is sufficient to acquire a vacuum 
image with the same imaging parameters as for the im- 
ages required for the reconstruction. This vacuum image 
can be used to normalize the three images of the specimen 
with respect to the illumination by a plane wave of ampli- 
tude 'one'. Subsequently the differences D^ ^/ 2 )( u ) be- 
tween the Fourier transformed 'conventional' image and 
each of the Fourier-transformed phase contrast images 
are calculated. Using Eq. (|3| and Q these differences 
can be evaluated for all u 7^ 0: 



D( Q ,±n/2)(u) =Io{u) 



<(± 



r/2)(") 



=#(u)**(0)J5(u, 0)exp(-ix(u))(l ± i) 
+^(0)^* {-u)E(u,0)exp(ix(u))(l T i) 

(10) 

with the envelope damping function E{u" ,u') and the 
phase distortion function x( u )- Again all non-linear con- 
tributions vanish as they are independent of the phase 
shift of the undiffracted beam. Since -D(o,±7r/2) i u ) are 
complex-valued the above system of equations can be 
solved for the complex quantity 



D 



(0,+tt/2) + i-D(0,-ir/2) 



2(l + i) 



(11) 



Correction for aberrations and negligence of the absolute 
phase of the undiffracted beam (setting the phase of the 
undiffracted beam equal to zero) leads to ty(u) ■ |v?(0)| 
for all u 0. The modulus of the undiffracted beam 
|^(0)| arises from the integration of the absolute square 
of 9(u) ■ I* (0)| over the whole Fourier space except for 
u=0: 



5>|*(0)| 2 . / ^{u'tfdu 1 



(12) 



and the intensity 1(0) of the Fourier transformed images 
at u=0. 1(0) can be evaluated using equation Eq. Q: 

7(0) = f \y(u')fdu' = |*(0)| 2 + f \^{u')\ 2 du'. 

(13) 

The above system of equations can be solved for |^(0)|: 



*(o)i = J km ± y/nor-^) (14) 



Since Eq. (12 1 and (13 1 are quadratic equations, two so- 



lutions exist, of which we found the negative signed to 



be the appropriate one. Using the quantity vf'(u) ■ |\&(0)|, 
obtained by Eq. (Ill, this finally leads to the Fourier 
transformed object wave function (u) for all u includ- 
ing u=0. The object wave function in real space \I/(r) 
is obtained by inverse Fourier transformation. Although 
the reconstruction of the object wave function in Fourier 
space is based on the same principle as the reconstruction 
in real space, it is quite sensitive towards the damping 
envelope functions at higher spatial frequencies, due to 
the aberration correction during the reconstruction pro- 
cess [2U]. An aperture function should be set to prevent 
a strong influence of the damping envelope at higher spa- 
tial frequencies. It will reduce the spatial resolution of 
the reconstruction slightly, because it has to be slightly 
smaller than the corresponding aperture of the spatial 
resolution of the original images. 



VERIFICATION OF THE RECONSTRUCTION 
TECHNIQUE 

To obtain a real test wave function for the reconstruc- 
tion procedure off-axis holography was used. The re- 
construction was carried out taking into account spatial 
frequencies up to 3 nmT 1 . The test sample consisted 
of nanoscaled Pt-particles on an amorphous (a-)C film. 
No reflections were visible in the sidebands of the holo- 
grams. The Pt-particles have a diameter of 1.5 — 2 nm 
which corresponds to a phase shift of ps 0.127T. The a- 
C carrier film induces a phase shift of about 2/37T with 
respect to the vacuum wave. The amplitude is 1 in the 
vacuum and decreases to about 0.85 for the carbon film. 
Based on these observations the Pt-particles can be con- 
sidered as weak-amplitude strong-phase objects. Strong 
reflections are not observed in the diffraction pattern, 
and the intensity of undiffracted electrons is much larger 
than any diffracted intensity. Linear image formation 
therefore applies to this particular test object. The ob- 
ject wave function obtained by off-axis holography will 
be denoted in the following as original wave function. 
Figure [2] shows three images which are calculated on the 
basis of the amplitude- and phase-distribution of the test 
object. The images were simulated assuming linear im- 
age formation and different phase shifts <f>pp = (Fig. [2] 
a)), 4> PP = tt/2 (Fig. gb)) and PP = -tt/2 (Fig. |]c)) 
in the wave aberration function. An ideal phase shifting 
device without a blocking structure was assumed. The 
images are calculated for a 200 keV transmission electron 
microscope with Z — —150 nm and C5 = 1 mm. For 
better comparison of the contrast generated by the phase 
plate, the contrast in all three images is displayed with 
the same minimum and maximum gray scale value as well 
as linear scaling between them. This results in a dark im- 
age for the case of 4> PP = (Fig. [2] a)) because the (a-)C 
film generates only weak contrast at such a comparatively 
low defocus value. On the other hand, if the phase-plate 



FIG. 2. Simulated images calculated on the basis of an object wave function of platinum particles on a thin amorphous carbon 
film acquired with off-axis electron holography, a) conventional image 4>pp = 0, b) 4>pp = n/2, c) <j>pp — —tt/2. (Cs = 1 mm, 
defocus Z = -150 nm, 200 keV) 



phase shift is ±7r/2 the film generates strong bright con- 
trast (Fig. [2]b)) or moderately bright contrast (Fig. [2] 
c)). The vacuum region in the upper left corner shows a 
constant grey value in all three images, and can be clearly 
distinguished from the film in Fig. [2]b), c). In Fig{2]a) 
the Pt-particles show dark contrast vs. the film, which is 
hardly visible due to the chosen contrast scaling. They 
show a slight dark contrast with respect to the bright 
(a-)C film in Fig. [2]b) and distinctly bright contrast in 
comparison to the a darker (a-)C film in Fig. [2]c). Due 
to the weak-amplitude character of the Pt-particles, the 
contrast between the film and the particles depends on 
the phase shift induced by the specimen but also on the 
phase-plate phase shift. For arbitrary phase objects like 
this specimen, the optimal contrast can be found any- 
where between and ±7r/2, in contrast to weak-phase 
objects, where the optimum contrast is always observed 
at phase shifts of ±7r/2. To be more specific, the contrast 
of the Pt-particles with respect to the (a-)C film changes 
drastically, because the film only shows little amplitude 
contrast, but induces a strong phase shift on the elec- 
trons, whereas the phase shift of the Pt-particles them- 
selves is only weak compared to that of the film. This is 
more obvious from the original phase distribution of this 
particular wave function which is shown in Fig. [3] a) . 

Using the reconstruction method outlined in the previ- 
ous chapter, the phase- and amplitude-distributions can 
be reconstructed and compared with the original wave 
function. Fig. [3] shows the original (Fig. [3] a)) and re- 
constructed phase distribution (Fig. [3] b)). With the 
exception of a constant phase shift, which results from 
setting the phase of the undiffractcd electrons to zero 
during the reconstruction, the phase distribution is well 
reconstructed from the three images. Line scans across a 
Pt particle are shown in Fig. [3] c) for the original wave 
function in gray and for the reconstructed wave function 
in black. Apart from the slightly lower resolution due to 
an aperture function, the phases are in very good agree- 
ment. 



For the verification of the reconstruction from images 
containing non-linear image contributions, another test 
object wave function is chosen. The object wave func- 
tion for single crystalline silicon along the [100]-zonc axis 
with a thickness of 10 unit cells corresponding to 5.431 
nm was calculated using the STEMsim program [2] . Va- 
cancies were introduced to generate differences in the 
local phase shift. The vacancies in the atom columns 
were introduced by removing silicon atoms in the phase 
gratings for the multislice algorithm. To simulate an oc- 
cupancy of 0.5 in a particular column, phase gratings 
of a perfect and defective Si(100) plane were used al- 
ternately. The super-cell contains in addition a vacuum 
region which is required for the reconstruction method 
in real space. Figure 4 shows images for different phase 
shifts of a) <fipp = 0, b) <j)pp = ir/2 and c) <j)pp = — tt/2, 
calculated using an equally spaced numerical integration 
algorithm tailored for phase-contrast imaging. A 200 keV 
transmission electron microscope with Cs = 0.5 mm and 
Z = 5 mm was assumed. Based on the three images 
the image wave function was reconstructed. The object 
wave function was retrieved by WTF-correction. Figure 
[5]shows the original (Fig. [5]a)) and reconstructed phase- 
distribution (Fig. [5]b)), as well as line scans (Fig. [5]c)) 
across an atom column containing vacancies and two ad- 
jacent fully occupied atom columns. The atom columns 
are located at regions with a large phase shift displayed 
'in white'. Apart from little noise which was simulated 
into the images (see next chapter), the reconstructed 
phase is in good agreement with the original phase used 
for the image simulation. It should also be noted, that 
due to the strong amplitude contrast in these images, 
the contrast between fully occupied atom columns and 
those with vacancies is barely visible by eye, but is re- 
constructed correctly, as shown in the line scans in Fig. 
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FIG. 3. : a) Phase of the original wave function acquired with off-axis electron holography scale from— -k (white) to +7r (black), 
b) reconstructed phase from the three images simulated in figure [2D . c) Line scans across a particle, taken from a) (gray line) 
and taken from b) (black line). The position is marked in a) and b) by the white dashed line. 
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FIG. 4. Simulated high-resolution images of silicon crystal 
along the [100]-zone axis containing vacancies (upper left cor- 
ner) calculated with TCC formalism for a) conventional image 
cj>PP — 0, b) <f>pp = 7r/2 , c) image with 4>pp — —tv/2 assum- 
ing a 200keV electron microscope with Cs = 0.5 mm and 
Z = 5nm 




FIG. 5. : a) Phase of the original wave function of the silicon 
sample used to simulate the three images in Figure |4j b) from 
the three images reconstructed phase of the wave function. 
White arrows point at the vacancy position, c) Line scans 
across two atom positions. Black taken from a) across the 
vacancy position and two adjacent atom positions, and gray 
taken from b) at the same location. Apart from little noise 
the line scans are in very good agreement, and show that 
the reconstruction works even if non-linear contributions are 
present in the simulated images. 



LIMITATIONS 

Due to the complexity of the equations used for recon- 
struction, one might expect that noise completely pre- 
vents the reconstruction of the object exit wave function. 
Therefore a series of image simulations was carried out 
to simulate the effect of different noise levels. The noise 
was simulated as Poisson noise for different average elec- 
tron doses per image pixel. The examination shows that 
only at lower electron doses below 500 electrons per pixel, 
a significant effect is observable. However, even at 250 
electrons per pixel the phase distribution in a crystalline 
silicon specimen can be reconstructed properly. Fig. [6] 
shows calculated peak signal-to-noise ratios (PSNR) from 
the original phase and reconstructed phase for different 
electron doses. The plot shows a drop in the PSNR at 
lower electron doses, approximately starting at 500 e/px. 
At this dose the phase signal starts to degrade, but is still 
in good agreement with the original phase down to noise 
levels where the SNR of the images taken is too small 
at all. While 20dB PSNR are considered acceptable for 
noisy images, values below this mark show typical arte- 
facts at high noise levels, which cause local deviations 
from the original distribution. Fig. [7] shows a series of 
reconstructed phases for Si in [100]-zone axis for differ- 
ent doses and the original phase for comparison. The 
top image with an electron dose as low as 250 e/px still 
shows a phase distribution which is in good agreement 
with the original phase, apart from being much noisier. 
250 electrons per pixel approximately translate to a SNR 
of 4 for the object wave function under investigation. To 
achieve a high SNR coherence is important. For strong- 
amplitude specimens, like almost all crystalline objects 
in zone-axis orientation, the dominant part in the exper- 
imental images is likely to be |^o| and and only 
little contrast remains for the phase part. It can be con- 
cluded that a high coherence should reduce the noise in 
the phase contrast part, which is necessary because it is 
the only information used for reconstruction. Non-linear 
terms add to this problem, therefore it is also important 
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1 °0 250(4.1:1) 500(5.6:1) 750(7:1) 
Dose [e/px] (SNR) 

FIG. 6. : Peak Signal-to-Noise Ratio (PSNR) calculated from 
reconstructed phase and original phase of a crystalline silicon 
sample as a function of the different electron doses per pixel. 
At low doses, significantly lower than 500 e/px, the phase 
signal degrades. Signal-to-Noise for 250 e/px, 500 e/px and 
750 e/px for the Si crystal specimen are given in the brackets. 
PSNR values above 20dB are acceptable for noisy images. 

that the dynamic range and number of grey levels of the 
camera (i.e. bit depth) is large enough for the specimen 
under investigation. For weak-phase objects this is not 
expected to be a problem but this behaviour was often 
observed for object wave functions with a strong am- 
plitude contribution if an 8-bit camera was assumed in 
the simulations with a limited number of discrete counts. 
Of course this depends significantly on aberrations and 
the spatial-frequency spectrum of the object wave func- 
tion. Under these conditions a correct reconstruction is 
impossible and can easily be identified by a very noisy 
reconstructed phase of the wave function. Coherence 
and SNR is expected to be one of the most important 
experimental problems of this method. Very small illu- 
mination angles for high spatial coherence and an un- 
diffracted beam with a small diameter requires a field- 
emission source electron microscope. A high mechanical 
and electrical phase-plate stability and specimen stabil- 
ity are needed because long exposure times are required 
for a high SNR. Another limitation arises due to the fact 
that real physical phase plates have different character- 
istics than assumed here for the reconstruction method. 
Depending on the type of phase plate, they introduce dif- 
ferent perturbations in the back focal plane of the objec- 
tive lens. In the case of a Boersch phase-plate, the elec- 
trode ring blocks a certain range of spatial frequencies. 




FIG. 7. : Reconstructed phase with different electron doses 
per pixel, a) 250 e/px b) 1250 e/px c) 2500 e/px d) original 
phase of the wave function. The comparison shows the degra- 
dation of the phase signal. As the peak signal-to-noise ratio 
calculations suggest, noise does not inhibit the possibility to 
reconstruct the wave function with the proposed method more 
than is expected. 



The electrode structure also induces a cut-on frequency 
at which diffracted electrons start to be phase shifted 
with respect to the undiffracted electrons. Other phase- 
plate types may have other characteristics depending on 
how they shift the phase of selected electrons. Simula- 
tions were performed with different phase-plate geome- 
tries and different shapes of electrostatic potentials in 
the back focal plane for different object wave functions. 
Their effect was taken into account by a special shape 
of the phase-plate function which is a complex function 
of u: P(u) = A(u)exp(i(j>pp(u)). It was found that the 
accuracy of the reconstruction depends significantly on 
how well the perturbation is confined to small spatial 
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frequencies u (near the undiffracted electrons) compared 
to the relevant information in Fourier space. The re- 
construction fails, if a significant fraction of electrons is 
diffracted into u values which are strongly perturbed by 
the physical phase plate itself. On the other hand the 
reconstruction works perfectly, if the diffracted informa- 
tion is contained mainly at spatial frequencies which are 
not affected by the physical phase plate. A good example 
for this behavior is a crystalline silicon specimen, where 
even with vacancies most of the information is contained 
in relatively large spatial frequencies which are not influ- 
enced by a perturbation in the vicinity of the undiffracted 
beam. This demonstrates that the influence of real phys- 
ical phase-plates can hardly be generalized. It always 
depends on the type of phase-plate and the size of the 
feature of interest. Electrostatic phase plates offer an sig- 
nificant advantage compared to carbon film phase plates 
if objects are studied which are not weak-phase objects. 
Variable phase shifts can be imposed, which is necessary 
for object wave reconstruction suggested in this work. 
For arbitrary wave functions, the reconstruction of the 
amplitude and phase requires three images with different 
phase shifts. Only two images can be obtained with one 
single carbon film phase plate, which is not sufficient for 
a proper reconstruction as demonstrated in chapter 3. It 
can also be understood from writing the complex function 
in sine and cosine functions which are not injective func- 
tions, and therefore the phase cannot be reconstructed 
uniquely from only two images. It is often discussed to 
what extent phase-plates increase information that can 
be obtained from a single image, mainly by applying 
the weak-phase-object approximation. This approach is 
commonly used due to its accessibility and experimen- 
tal importance to justify experimental efforts made in 
phase-plate development. Recently, higher order series- 
expansion has been applied to investigate strong-phasc- 
objects and the increase in contrast of the obtained phase 
contrast images [ST]. Carbon-film phase-plates with dif- 
ferent thicknesses are suggested. However, for arbitrary 
wave functions, which are of concern to material science, 
one will immediately conclude, that using a phase plate 
with a fixed phase shift, whatever phase shift that would 
be in particular, will result in strong phase contrast or 
none at all - very much depending on the object wave 
function itself. 



SUMMARY 

Phase-shifting devices with the capability of imposing 
variable phase shifts can be used to reconstruct the am- 
plitude and phase of wave functions from arbitrary ob- 
jects on the basis of only three different images. Our 
method relies on the difference between two images ac- 
quired under the same experimental conditions apart 
from the phase shift. This effectively eliminates all con- 



tributions due to non-linear image formation which can- 
not be avoided for most materials science specimens un- 
der well-defined diffraction conditions. Complex equa- 
tions or methods are not required to account for non- 
linear image formation. For a transmission electron mi- 
croscope without corrector, the object wave function can 
be as usual obtained from the image wave function by 
taking into account the wave transfer function. A re- 
liable reconstruction depends on the signal-to-noise ra- 
tio. Using silicon as a test object it was shown that 
a proper reconstruction is possible for electron doses as 
low as 250 electrons/pixel. The resolution of the recon- 
structed wave function corresponds to the resolution of 
the images, which is only marginally reduced during the 
reconstruction process. 
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